A phase-separation perspective on dynamic heterogeneities in glass-forming liquids 
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We study dynamic heterogeneities in a model glass-former whose overlap with a reference config- 
uration is constrained to a fixed value. The system phase-separates into regions of small and large 
overlap, so that dynamical correlations remain strong even for asymptotic times. We calculate an 
appropriate thermodynamic potential and find evidence of a Maxwell's construction consistent with 
a spinodal decomposition of two phases. Our results suggest that dynamic heterogeneities are the 
expression of an ephemeral phase-separating regime ruled by a finite surface tension. 



The conspicuous lack of a growing correlation length, 
contrasting with the very steep increase of the relaxation 
time, has been a puzzle in the physics of structural glasses 
for quite a long time. Arguably, the first breakthrough 
has been the discovery of dynamic heterogeneities fH , and 
the detection of a growing dynamical correlation length, 
[3, S] ■ If we take two snapshots of the system sepa- 
rated by a time lag comparable to the a relaxation time, 
Tq., the particle displacements vary enormously across the 
system, and the typical size of the mobility-correlated 
regions increases on lowering the temperature. 

More recently, by studying the thermodynamics of sys- 
tems subject to amorphous boundary conditions an 
entirely different, fully static, correlation length has 
been discovered [1, 0] ■ Cs also grows upon cooling, even 
though its surge occurs at lower temperatures than ^d- 
The static correlation length has a natural interpreta- 
tion as the size of the cooperatively rearranging regions 
Q , and within the random first-order theory Q it is de- 
termined by the balance between a surface tension cost 
and a configurational entropy gain of a rearrangement. 

Clearly, it would be desirable to unify the dynamic and 
the thermodynamic frameworks, so as to understand the 
interplay between the two correlation lengths. Although 
the static-dynamic connection is clear in mean-field sys- 
tems [iq| and so me p rogresses have been made in more 
realistic systems ll|, we are quite far from a unifying 
picture in real glass- formers. Here we show that surface 
tension, which is a crucial ingredient of the thermody- 
namic framework, also plays a key role in the formation 



of dynamic heterogeneities. In so doing, we establish a 
much-needed further link between dynamic and thermo- 
dynamic relaxation in glass-forming liquids. 

Let us start with the standard measurement of the 
dynamic correlation length ^d- Our glass-former is the 



well-known soft-sphere model in S-d [26|. A useful tool 
to measure ^d is the overlap, which quantifies how much 
a configuration at time t is similar to the reference con- 
figuration at t = 0. If we partition the system in small 
cubic boxes and let ni be the number of particles in box i, 
the local overlap is defined as q(ri,t) = ni{t) ni{0), where 
Ti refers to the centre of cell i [27|. The spatial map of 
the local overlap tells us how much different regions of 
the system have decorrelated (respect to the initial con- 
figuration) over a time t. In Fig. [1] (top) we show two 
snapshots of the overlap field. We see that a,t t = Ta 
there are large heterogeneous regions, which eventually 
fade away for longer times. To quantify their size we 
must compute the overlap correlation function, 

G(r,t) = ((7(0,t)<7(r,<)) - {q{0,t)){q{r,t)), (1) 

or its Fourier transform S{k,t) (Fig. [2l left). In general, 
given a correlation function in Fourier space, it is well- 
established practice to extract the correlation length 
^ from the small-fc linear interpolation of S^^ vs. fc^, 



S{k,ty 



A + Bk^ 



(2) 
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from which the correlation length is obtained as £,{1)^ — 
B/A. The validity of Eq.[2]is shown in the inset of Fig. 2, 
right. ^28i]. The time-dependent correlation lenght £^{t) 
represents the size of the dynamical heterogeneities at 
time t. This lengthscale grows as the time approaches 
Ta and the heterogeneities become more extended (inset 
of Fig. 2, left). The largest value of £,{t) (reached at the 
Ta) defines the so-called dynamical correlation length. 
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FIG. 1: Fluctuations of the overlap field, 5q{r,t) = q{r,t) — 
{q{t)) for a 2-d slice of the system. Upper panels: uncon- 
strained system. Lower panels: constrained system (Q — 
0.25). Left panels: t — Ta. Right panels: large times. L — 16. 
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FIG. 2: S{k, t) at different times for the unconstrained (left) 
and constrained (right, (Q = 0.25) cases. Left inset: cor- 
relation length ^(t) as extracted from Eq. Right inset: 
S{k,t)-'^ vs. k^. T = Tmc and L = 16. 



What happens beyond r^? The memory of the initial 
configuration is gradually lost, so that the correlation 
function S{k,t) decays sharply, Fig. 2 (left). What hap- 
pens to £,{t) is less clear, because the vanishing of S{k,t) 
makes it hard to fit a reliable value of ^{t) through Eq. [2] 
Although this point is debated 14], our results indicate 
that ^{t) decreases beyond Tq, (left inset of Fig. 2), in 
line with other studies [13]. However, what is important, 
and definitely out of question, is that the correlation de- 
creases for large times, irrespective of its spatial range. 
Heterogeneities blur as q{r,t) becomes zero everywhere, 
and their size ^(t) becomes somewhat ill-defined. 

Let us now make a different experiment. We want to 
impose a constraint on the dynamics, so that the sys- 
tem cannot entirely lose memory of its initial configu- 
ration. This can be implemented by imposing a lower 
bound on the global overlap, Q{t) = 1/V J dr q{r,t). In 
the unconstrained case Q{t) goes asymptotically to zero, 
as the memory of the initial configuration fades [29*] . On 
the other hand, if we run the dynamics with the con- 
straint Q(t) > Q things change [30|. Initially the system 
does not feel the constraint: the global overlap decreases 
from its t = value, Q — 1, and everything proceeds 
as described above, including the growth of the hetero- 
geneities. However, at later times Q{t) hits its lower 
bound Q and it cannot decrease further. What happens 
to the dynamical heterogeneities in this case? 

There are two alternative hypotheses. First, the corre- 
lation S{k,t) and its spatial range ^{t) decay to zero for 
large t as in the free case. Due to the constraint, however, 
such endgame cannot happen in the same manner as in 
the free case, i.e. with q{r,t) becoming zero everywhere. 
Heterogeneities must thus become very small, forming a 
salt-and-pepper configuration of the field q{r,t), so that 
the total integral of the field stays equal to Q. Yet, 



if there is a nonzero surface tension between high and 
low overlap regions, such a scenario is not what we ex- 
pect: the surface tension would force different domains 
to merge, driving the system towards a phase-separated, 
highly correlated state 15|. Hence, the second hypoth- 
esis is that the correlation does not decay and that the 
dynamic correlation length ^(i) grows beyond ^d, up to 
an asymptotic value of the order of the system's size L. 
The stark difference between these two hypotheses sug- 
gests that the constrained experiment may clarify the 
mechanisms of formation of dynamical heterogeneities. 

Inspection of the overlap field in the constrained case. 
Fig. [1] (bottom), is quite telling: for large t the system 
phase separates into high and low overlap regions, form- 
ing stable dynamical heterogeneities of the order of the 
system size. From a quantitative point of view the sit- 
uation is equally clear: in contrast with the free case, 
the constrained correlation function does not go to zero 
for large times, but saturates at a finite value. Fig. [2] 
(right). Hence, even in the late time regime dynamic 
heterogeneities remain strongly correlated. 

The study of the correlation length in the constrained 
case confirms this scenario [3l|. In an infinite system 
undergoing phase separation, or at a critical point, the 
intercept A in Eq. [2] vanishes while the slope B remains 
finite, so that ^ grows indefinitely 1^. On the other 
hand, in a finite system phase separation means that 
^ becomes comparable with system size L. As a con- 
sequence, the finite-size (periodic) real space correlation 
function G{r, t) ceases to be a simple exponential for large 
r (small k). This implies that the intercept A in Eq.[2]can 
go below zero and take small negative values (0(1/L^)). 
In a finite-size system it is therefore convenient to com- 
pare A/B = vs. to check whether or not phase 
separation occurs. From Fig. [3] (left) we see that in the 
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FIG. 3: Left: A/B = C'^^t) (see Eq. EJ at T = Tmc in 
the constrained (circles) and unconstrained (triangle) cases. 
Right: the same at T = I.SSTmc- L = 16. 
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FIG. 4: Energy difference AE{t) = E{t)~Eo vs f at T = Tmc 
with constrained dynamics, Q = 0.25. Eq is a parameter of 
the fit E{t) = Eo + 7t"^/^. The line is corresponding 
to the surface tension exponent 9 = 2. 



unconstrained case keeps well clear of L^^, while 
in the constrained case £,~'^{t) unmistakably goes below 
L~^. This is exactly what we expect in a system with 
nonzero surface tension undergoing phase separation. We 
studied two other sizes, L — 8 and L — 25, and in both 
cases (,~'^{t) drops below L^^, indicating phase separa- 
tion. At higher temperatures, however, though the cor- 
relation is enhanced by the constraint, the latter is in- 
effective to make ^~^(i) drop below (Fig. |31 right). 
These results are consistent with the idea that the sur- 
face tension decays at high temperature, thus preventing 
phase separation [l6l |. 

In systems with conserved order parameter undergo- 
ing phase separation the domains size ^(i) grows as t^/"^ 
and the dynamics proceeds by reducing the total amount 
of interfaces, and therefore of energy, in the system [l5| . 
The interface energy per domain scales like where 9 
is the surface tension exponent. The total number of do- 
mains is L'^/£^'^, so that the total interface energy density 
is AEit) - ^j^g standard case 

6* = d - 1, so that AE{t) - 1/^1/3 ^ Figg] shows 
that something remarkably similar happens in our case. 
After the constraint kicks in, AE{t) decays compatibly 
with an exponent 1/3. Hence, even though fitting coars- 
ening exponents is notoriously difficult, and one must be 
careful in drawing any conclusion, our data seems to be 
compatible with the 'naive' exponent 9 = 2 3 121- 

In general, phase-separation is the landmark of first 
order phase transitions and metastability. At the mean- 
field level one can normally define a thermodynamic po- 
tential as a function of the order parameter that, below 
some spinodal point, exhibits a stable and a metastable 
minimum, corresponding to the two phases. In finite di- 
mension Maxwell's construction makes the potential con- 
vex, so that the derivative of the potential is constant 
(zero second derivative) in a finite interval (Fig. 5, in- 
set). Maxwell's construction implies that when the order 



parameter is conserved and constrained to take a value 
in the non-convex interval, phase separation occurs. We 
have clearly observed phase-separation. Can we define a 
thermodynamic potential displaying Maxwell's construc- 
tion? 

Let us proceed minimalistically. Our phase-separating 
order parameter is the overlap Q, so it is a potential 
W{Q) we are after. Besides, the potential must deter- 
mine the observed probability distribution of Q through 
the relation, P(Q) = exp[-NW{Q)] 9(0 - Q), where 
the 0-function enforces the constraint [18|. If we com- 
pute the average linear fluctuation of Q and expand the 
exponential, we obtain. 



(3) 



This quantity is easy to compute: we let the system 
evolve until the constraint is hit, and then we measure 
the (very small) average fluctuation of the overlap Q 
over Q (see [l^l for a different definition of the poten- 
tial). We report W'{Q) in Fig. ([5]). The second deriva- 
tive of the potential is clearly nonzero at high T , whereas 
around the Mode Coupling temperature a finite region 
with W"{Q) ^ develops. This is evidence of Maxwell's 
construction and it supports the link between phase sep- 
aration and metastability in our system. 

W{Q) is a finite-dimensional variant of the two-replica 
potential originally introduced in mean-field spin-glasses 
[20| , and later generalized to structural glasses [2l| . This 
potential is the free energy cost to keep a configuration 
(the running one in the present work) at fixed overlap 
Q with a generic equilibrium configuration (the initial 
reference one). Below a dynamic transition (roughly, the 
Mode Coupling temperature), the mean- field potential 
develops a metastable minimum at a finite value of Q. 
In this framework relaxation at low temperatures can 
be interpreted as a barrier crossing process, bringing the 
system from the metastable minimum (short times, finite 
Q) to the stable minimum (long times, zero Q) [22]. The 
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FIG. 5: The derivative W'{Q) of the thermodynamic po- 
tential at different temperatures, ranging from 2.13Tmc to 
0.89rjv/c- Inset: a cartoon of Maxwell's construction for the 
potential and its derivative. 



constraint is just a stratagem to keep the overlap within 
the nonconvex region of the potential, as to interrupt 
relaxation and therefore force phase separation. 

We have studied dynamic heterogeneities in a glass- 
forming liquid with constrained global overlap. At low 
temperature both the dynamic correlation function and 
the thermodynamic potential indicate that the system 
phase separates into regions of high and low overlap. On 
the contrary, at high temperature no phase separation 
occurs, supporting the view of a surface tension that de- 
creases at high T. The co-existence of regions belonging 
to different amorphous 'states' (here the high/low overlap 
patches) is reminiscent of the random first-order theory 
of thermodynamic relaxation [9] . In the dynamical case 
the evolution of these regions is driven by a classic coars- 
ening mechanism, which is stable with the constraint, but 
ephemeral in absence of the constraint. In the thermo- 
dynamic case, on the other hand, the evolution of these 
regions is presumably driven by an entropic mechanism 
0. Our results show that surface tension and metasta- 
bility stand as key links between the two frameworks. 
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